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Abstract 

The energy dependence and intensity of Coulomb interaction between quasiparticles in metallic wires is obtained 
from two different methods : determination of the temperature dependence of the phase coherence time from the 
magnetoresistance, and measurements of the energy distribution function in out-of-equilibrium situations. In both 
types of experiment, the energy dependence of the Coulomb interaction is found to be in excellent agreement with 
theoretical predictions. In contrast, the intensity of the interaction agrees closely with theory only with the first 
method, whereas an important discrepancy is found using the second one. Different explanations are proposed, 
and results of a test experiment are presented. 
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1. Introduction 

The description of electrical transport in metals 
is based on the existence of long-lived quasipar- 
ticles. The finite quasiparticlc lifetime appears in 
mesoscopic physics as a limitation of their phase 
coherence time, which determines the amplitude 
of quantum interference effects. The three kinds 
of processes that limit the quasiparticle lifetime 
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in metals are electron-phonon scattering, electron- 
electron scattering [1], and spin-flip scattering of 
electrons from magnetic impurities [2,3]. At tem- 
peratures below about 1 K, the rate of electron- 
phonon scattering is weak, and in metallic samples 
without magnetic impurities the dominant inelas- 
tic scattering process should be the Coulomb in- 
teraction between electrons [1] . 

In this paper, we focus on experiments per- 
formed on very clean (99.9999%) silver wires, in 
which the effect of magnetic impurities is expected 
to be small [4,5]. We review the results obtained 
from weak localization measurements, in which 
the phase coherence time t v (T) is extracted, and 
from energy relaxation experiments, in which the 
energy exchange rate between quasiparticles is de- 
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rived from their energy distribution function f(E). 
In the former experiments, we find that both the 
temperature dependence and overall magnitude of 
t v (T) agree with the theoretical predictions. In the 
latter experiments, the energy dependence of the 
inelastic rate agrees with theoretical predictions, 
but the overall magnitude fluctuates significantly 
from sample to sample. 



2. Two experiments for measuring 
Coulomb interaction between QPs 

In metallic thin films, quasiparticles (QPs) expe- 
rience frequent elastic scattering from grain bound- 
aries, film edges and impurities. In this diffusive 
regime, characterized by a diffusion constant D, 
the screening of the Coulomb interaction is re- 
tarded, and the corresponding (squared) matrix el- 
ement between two QPs, derived by Altshuler and 
Aronov in the early 80s [1], depends on the en- 
ergy e exchanged during the interaction process: 
|M(e)| 2 oc e~ 3 / 2 in quasi-one-dimensional wires. 
This energy dependence results in a temperature 
dependence of the phase coherence time t v (T) oc 
jgj^ w } uc ] 1 h as been observed in aluminum 
and silver wires by Wind et al. [7] down to IK, 
and by Echternach et al. [8] in gold wires down 
to lOOmK. The most convenient method to access 
t v is the measurement of the magnetoresistance of 
wires with a length L long compared to the phase 
coherence length L v = -J ' Dt Vi which exhibits a 
small peak or dip at zero magnetic field due to weak 
localization [9]. When the rate of spin precession 
due to spin-orbit coupling exceeds the dephasing 
rate, as is usually the case at low temperature, the 
relative amplitude of the zero-field dip in the resis- 
tance gives direct access to L v : 

SR ^ 2R L v 

with R K — h/e 2 rts 26 kfl the resistance quantum. 
The width in field of this dip corresponds to a flux 
quantum in the area L v w, with w the wire width. 
In practice, magnetoresistance curves measured at 
different temperatures are fit with a theoretical ex- 
pression for ^jf{B) in which the only fit parameters 
are the phase coherence length L v , the spin-orbit 



length L so , and the width of the wire w [4]. The 
two last parameters, L so and w, are fixed at a con- 
stant value independent of temperature for each 
sample [10]. Then, t v is obtained as L 2 /D, with D 
obtained from the resistance R = — K-n-^j where 
vp is the density of states at the Fermi energy (2 
spin directions) and t the wire thickness. In order 
to compare with theory, the resulting curve t v (T) 
is fit with 

t v (T) = (AT 2 / 3 + BT 3 )- 1 . (1) 

where AT 2 / 3 is the Coulomb interaction rate and 
BT 3 the approximate electron-phonon scattering 
rate [11]. 

In theory, the exchange part of the Coulomb in- 
teraction leads to [12] 

A = ±(^£-*X'* (2) 
h \Av F LwtR K J y ' 

The contribution due to the Hartree term has not 
been evaluated for wires [13]. 

Another experimental method to access the in- 
teraction processes consists in driving the QPs out- 
of-equilibrium by a finite voltage U between two 
contacts at the ends of the wire, which act as QP 
reservoirs [14]. At energies between — ell and 0, 
the diffusion of QPs from the occupied states at 
one end to empty states at the other end results, 
in absence of inelastic processes, in a two-step dis- 
tribution function f x (E) inside the wire as pic- 
tured in Fig. 1. (The shorthand f x (E) stands for 
f(x, E), where we measure distance in units of the 
wire length L, so that < x < 1.) This distri- 
bution function can be understood as a linear in- 
terpolation between the distribution functions at 
the boundaries of the wire. Electron-electron inter- 
actions lead to a redistribution of energy between 
QPs at each position, hence to a rounding of f x (E) . 
In experiments, f x (E) at a given position in the 
wire is deduced from the differential conductance 
dI/dV(V) of a tunnel junction between a super- 
conducting probe electrode and the wire. In order 
to relate f x (E) to the matrix element of the inter- 
action, the data are fit with the solution of the sta- 
tionary Boltzmann equation in the diffusive regime 
[15,16]: 
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f x (E)n x (E) 




Fig. 1. Schematic diagram showing the spatial and en- 
ergy dependence of the distribution function f x (E) of QPs 
driven out-of-equilibrium by the voltage U using the ge- 
ometry of Fig. 6 with the switch in position 1. The sur- 
rounding box shows the uniform density of states in the 
metal and the gray volume shows the occupied states whose 
normalized density is f x {E)n x (E). The thick line shows a 
typical double step distribution function at x = 1/4 as in 
Fig. 7. 

1 d 2 f x (E) 
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-T* xl (x,E,{f})-T^ x (x,E,{f})=Q 

(3) 

where 2^ (x, E, {/}) and 1°^ (x, E, {/}) are the 
rates at which quasiparticles are scattered in and 
out of a state at energy E by inelastic processes. 
The diffusion time tb = L 2 /D is the typical time 
spent by a QP in the wire. Assuming that the dom- 
inant inelastic process is Coulomb interaction be- 
tween QPs and phonon emission or absorption, the 
inelastic scattering integrals read 

Z^S (x, E, {/}) = J de f x (E) (1 - UE - e)) W(e) 

X^a (x,E,{f}) = J de f x (E + e)(l-f x (E))W(e) 

with 

W(e) = W e ^ e (e)+W e . ph (e) 

W e - e (e) = K(e) j dE'f x (E')(l - f x (E' + e)) 

W e - ph (e) = n ph s 2 (n ph (\e\) + 0(e)). 

The kernel function K (e) = n ee e~ 3 ^ 2 is pro- 
portional to the averaged squared interaction ma- 
trix element |M(e)| 2 between two quasiparticles 
exchanging an energy e [1]. Its intensity n ee , which 
can be derived either from the expression of the 
microscopic interaction potential [17,18], or from 
the fluctuation-dissipation theorem [18], is [19] 



K ee = (y2DirH 3/2 v F wtj . (4) 

This derivation takes into account the exchange 
term only. The Hartree contribution to K (e) is 
expected to be smaller [1,17]. The electron-phonon 
coupling has an intensity n p h and is proportional 
to the sum of the Bose energy distribution of 
phonons n p /j(|£|) representing stimulated absorp- 
tion or emission of phonons and the Heaviside 
function 9(e) representing spontaneous emission. 
A more accurate description of electron-phonon 
coupling was developed in [11]. However, we re- 
strict here to the simplistic form for W e - P h be- 
cause the effect of phonons is very small. Thus, 
for all the fits to the experiments, we fix the value 
of K p h at 4 ns ^^meV 3 , which is compatible with 
the weak localization measurements [2 2]. 

The boundary conditions for Eq. (3) are Fermi- 
Dirac distributions at the ends of the wire, with a 
temperature higher than the cryostat temperature 
due to electron heating in the reservoirs [23,20,24]. 

The link between the two parameters determin- 
ing the effect of Coulomb interaction, A and n ee , 
can be made explicit by noting that the dephasing 
rate is the average of the inverse of the lifetime of 
QPs at energies within fcgT of the Fermi energy 
[25]: 
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so that 
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4K ee k B 
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While this derivation reproduces the correct de- 
pendence on sample parameters of the more rigor- 
ous theory [6,12], the prefactor depends on the ex- 
act value of the cutoff, whose order of magnitude 
is hi Tip. The choice of the cutoff can be made so 
that our derivation stays consistent with the ex- 
pressions Eq. (2), (4) of A and K ee . Thus it is pos- 
sible to express A as an intensity ka for Coulomb 
interaction, using 
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3. Comparison between experimental and 
theoretical results for both methods 

We present here data taken on wires deposited 
from 6N-purity (99.9999%) silver sources. The 
fabrication procedure for weak localization type 
(WL) samples is described in Ref. [4]. The sample 
parameters are given in Table 1 (weak localiza- 
tion measurements) and Table 3 (energy relax- 
ation measurements). The names of the samples 
used in energy relaxation (Relax) experiments 
contain Roman numerals, which indicate the in- 
dex of the experiment, and a number, which is the 
approximate wire length in microns. Most Relax 
samples were obtained in a single step, using two- 
angle evaporations through a suspended mask [20] . 
Samples AgII5 and AglllO, on the one hand, and 
AgIV20a and AgIV20/3, on the other hand, were 
fabricated at the same time, on the same chip. 
Samples AgXIlO, AgXII40 and AgXV40 were fab- 
ricated in two steps of e-beam lithography: in a 
first step, the wire pattern was defined, then sil- 
ver was evaporated and followed by a lift-off, and 
a new deposition of resist. In a second step, the 
pattern for the aluminum electrodes was exposed 
to the electron beam. In the vacuum chamber 
of the deposition machine, the silver layer was 
cleaned by argon ion milling. A thin (3 nm) layer 
of aluminum was then deposited, followed by an 
oxidation in 1.3 mbar of oxygen-argon (20%-80%) 
during 8 minutes, in order to form the tunnel bar- 
rier. Finally, a layer of aluminum was deposited. 

In Fig. 2, we present t v (T) for the first three 



Sample 


L w t R D 
((im) (nm) (nm) (kf2) (cm 2 /s) 


Ag(6N)a 
Ag(6N)b 
Ag(6N)c 
Ag(6N)d 


136 65 47 1.44 117 
271 100 45 3.30 69.2 
400 105 53.5 1.44 187 
285 90 36 2.00 167 



Table 1 

Geometrical and electrical characteristics of samples for 
weak localization measurements. The diffusion coefficient 
D is obtained using Einstein's relation 1/p = vpe D with 
the density of states in silver up = 1.03x 10 47 J _1 m~ 3 , and 
the resistivity p extracted from the resistance _R, thickness 
t, length L and width w of the wire. 




T(K) 

Fig. 2. Phase coherence time vs temperature in samples 
Ag(6N)a (■), Ag(6N)b (▼), and Ag(6N)c (•), all made of 
6N sources. Continuous lines are fits of the data to Eq. (1). 
The quantitative prediction of Eq. (2) for electron-electron 
interactions in sample Ag(6N)c is shown as a dashed line. 

WL samples (the data points of the last one, which 
are presented in Ref. [4], are so close to those of 
the third one that they would confuse the figure), 
as well as the best fits with Eq. (1). The fit pa- 
rameters are given in Table 2. The fit value of 
A is very close to the theoretical value for the ex- 
change contribution of the Coulomb interaction, as 
can be seen in Fig. 4 where the X-coordinate of 
the solid squares is the theoretical value of ka us- 
ing Eqs. (2) and (5), and the Y-coordinate is the 
value from experiment. 

The situation is quite different in energy relax- 



Sample 
■ 


A th y A B 

(ns-^- 2 / 3 ) (ns^K" 3 ) 


Ag(6N)a 
Ag(6N)b 
Ag(6N)c 
Ag(6N)d 


0.55 0.73 0.045 
0.51 0.59 0.05 
0.31 0.37 0.047 
0.47 0.56 0.044 



Table 2 

Theoretical predictions of Eq. (2) (A thy ) and fit parameters 
(A and B) for t v (T) in the samples of Table 1 using the 
functional form given by Eq. (1). Comparison of A thy and 
A is shown graphically in Fig. 4. 
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Sample 


L 


w 


t 


R 


D 


TD 




(jim) (nm) ( 


nm) 


(0) 


(cm 2 /s 


)(ns) 


AgI5 


5.05 


90 


43 


41 


121 


2.1 


AgII5 


5.2 


66 


39 


44 


173 


1.6 


AglllO 


10.3 


65 


39 


81 


191 


5.6 


AgIII20 


19.6 


160 


43 


45 


241 


16 


AgIV20« 


19.7 


95 


44 


86 


208 


19 


AgIV20/3 


19.9 


100 


44 


91 


188 


21 


AgX20 


21.7 


100 


48 


80 


214 


22 


AgXIlO 


9.55 


124 


45 


31 


211 


43 


AgXII40 


38 


180 


45 


108([26]) 


165 


87 


AgXV40 


38 


145 


45 


134 


165 


87 



Table 3 

Geometrical and electrical characteristics of samples for 
energy relaxation measurements. 




Sample 


^ee 




o 


(ns- 1 


meV- 1 ^) 


AgI5 


0.060 


0.95 


AgII5 


0.076 


0.5 


AglllO 


0.073 


0.54 


AgIII20 


0.024 


0.5 


AgIV20o 


0.043 


0.40 


AgIV20/3 


0.043 


0.37 


AgX20 


0.037 


0.11 


AgXIlO 


0.032 


< 0.18 


AgXII40 


0.025 


0.18 


AgXV40 


0.031 


0.32 



1 E/eU 



Table 4 

Theoretical predictions of Eq. (4) («;'e y ) and fit parame- 
ters («ee) for fx(E) in the samples of Table 3 using the 
solution of the Boltzmann equation Eq. (3). The distribu- 
tion functions measured on sample AgXIlO were so close 
to the nonintcracting regime that it was only possible to 
give an upper bound to the value of K ee . Comparison of 
4e y and K ee is shown graphically in Fig. 4. 

nearly an order of magnitude larger than the value 
given by Eq. (4). Similar discrepancies exist for the 
other Relax samples. It could be argued that the 
numerical prefactor in Eq. (4) is incorrect. Fig. 4 
seems to rule out this explanation: the circles cor- 
responding to the theoretical and fit values, given 
also in Table 4, present a large scatter, and so the 
ratio between experiment and theory does not ap- 
pear to be constant. 



Fig. 3. Measurements (o) and fits (solid curves) of the 
quasiparticle energy distribution function fi (E) for five 

2 

different values of the applied voltage U across the wire 
AgIV20a. The data have been shifted vertically for clarity. 

ation experiments. We show in Fig. 3 distribution 
functions f(E) measured in the middle of sam- 
ple AgIV20a, for U ranging from 0.1 to 0.5 mV, 
plotted as a function of the reduced energy E/eU. 
Solid lines are fits resulting from the numerical so- 
lution of the Boltzmann equation, obtained with 
K ee — 0.40 ns _1 meV -1 / 2 . The increase in slope of 
the middle step of f(E) when U increases, char- 
acteristic of the effect of Coulomb interaction, is 
well reproduced. However, the fit value for K ee is 



4. Discussion of the discrepancy between 
the two experiments 

Figure 4 reveals a very puzzling difference be- 
tween weak localization (WL) and energy relax- 
ation (Relax) experiments. Whereas the results of 
both types of experiments are precisely accounted 
for by the theory of Coulomb interactions in dis- 
ordered wires as far as the energy dependence is 
concerned, the prefactor is well understood for the 
first, but not at all for the second. In order to re- 
solve this puzzle, we now list the differences be- 
tween the two types of experiments: 
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theory (ns" 1 meV" 1/2 ) 

Fig. 4. Comparison of the experimental prcfactor with the 
theoretical prediction Eqs. (2), (4), (5), for weak localization 
experiments (■) and energy relaxation experiments (o). If 
we assume that a small amount of magnetic impurities is 
present in the WL samples, the fit values of Kj\ can be 
reduced down to the bottom of the dashed lines below the 
squares. Similarly for the Relax experiments, if we assume 
that magnetic impurities are present, we obtain a range of 
values of K ee compatible with the data, represented as a 
dashed line below the o. The behavior of sample AgX20 
was measured in a magnetic field, allowing us to place an 
upper bound on the concentration of magnetic impurities, 
and hence to rule out the possibility of reducing the value 
of Kee more than 15%. Thus, this data point is represented 
as a bold circle without any dashed line. 

4.1. Possibility of extrinsic energy exchange 
processes in Relax samples 

WL experiments are extremely sensitive to very 
small quantities of magnetic impurities. It was 
shown in [4] that even in our cleanest Ag(6N) 
wires, there was evidence for magnetic impuri- 
ties at concentrations of about 0.01 ppm, i.e. 1 
impurity atom for every 10 8 Ag atoms. Their 
contribution to t v was visible only at the lowest 
experimental temperatures. In Fig. 4, we have in- 
dicated with the vertical dashed lines how far the 
fit values of ka can be reduced if one includes a 
small concentration of magnetic impurities as an 
extra fit parameter. 

It is now undesrtood that magnetic impurities 



also mediate energy exchange between electrons 
[3,5]. Could the presence of magnetic impurities 
explain the anomalously large apparent values of 
n ee observed in many Relax experiments? Since 
most of the Ag samples used in the WL experi- 
ments were fabricated in the same deposition sys- 
tem used for the Relax samples, we expect that 
Relax samples should be equally clean. This hy- 
pothesis must be checked, however. The presence 
of magnetic impurities in Relax samples can be de- 
tected directly by performing the experiment as a 
function of magnetic field [5]. In samples AgX20 
and AgXIlO, the magnetic field dependence of the 
measurements set an upper bound to the concen- 
trations of magnetic impurities at 0.1 and 0.6 ppm 
respectively. For sample AgX20, if we include the 
effect of 0.1 ppm of magnetic impurities into the 
analysis of the Relax data, the value of 
duced by only 15%. In sample AgXIlO, the distri- 
bution functions were so close to the noninteract- 
ing regime that it was only possible to place an 
upper bound on K ee , hence this sample does not 
appear in Fig. 4. 

For the Relax samples that were not measured 
in a magnetic field, no upper bound to the concen- 
tration of magnetic impurities is experimentally 
determined. We have estimated the resulting sys- 
tematic uncertainty in n ee by the following anal- 
ysis. We have assumed that electron-electron in- 
teractions mediated by magnetic impurities con- 
tribute to energy exchange. For this process, the 
interaction kernel is approximately K(e) = K2£~ 2 
[3,27]. If we fit the data using the value of k 2 as 
an additional fit parameter, we can ask how small 
the value of become before the fits become 

clearly incompatible with the data. The results are 
shown by the dashed lines descending below the 
points for the Relax samples in Fig. 4. As can be 
seen, for some samples the fits are somewhat in- 
sensitive to the relative weights of K ee and K2, and 
the discrepancy between theory and experiment 
gets smaller. Nevertheless, the discrepancy still re- 
mains. We conclude for the time being that extrin- 
sic energy exchange processes with K{e) cx e~ 2 
are unlikely to explain completely the discrepancy 
between experiment and theory. This issue will be 
discussed further in section 6. 
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4.2. Sample dimensionality 

The intensity and energy dependence of 
Coulomb interaction depends on sample dimen- 
sionality [1]. The one-dimensional (ID) regime 
described in section 2 corresponds, in WL exper- 
iments, to situations where w,t <C L v <C L. This 
inequality is well obeyed in our experiments, where 
L v varies between 1 /im to 20 /im. In practice, 
the wire length L was chosen much larger than 
L v (r m i n ), where T m j n is the lowest experimental 
temperature, in order to reduce the amplitude of 
conductance fluctuations, which spoil the analysis 
of the magnetoresistance in terms of the WL the- 
ory. 

In Relax experiments, on the other hand, the dis- 
tribution function f(E) only contains information 
on the interaction process if it is far from a Fermi 
function and far from a perfect double-step, i.e. 
if L ss few L v (e/7 max /fcs). Thus the wire length 
is smaller than for the WL experiments. The di- 
mensionality criterion for Relax is illustrated in 
Fig. 5, where we plot the function K(e) calculated 
using the discrete sum over the longitudinal and 
transverse wave vectors [17,28] 
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271 
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m 
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K (e) oc V . 



(6) 



Qy,Qz 



where q x = q y = ^ and q z = ^ are 

the wave vector components with n x e N* and 
n y ,n z e N. 

Typical sample dimensions were chosen: L = 
10 /im, w = 130 nm, t = 45 nm and D = 
200 cm 2 /s. Fig. 5 shows that for all relevant en- 
ergies in the experiments, K{e) is far from the 
1D-3D transition. For small energies near /csT mi „, 
the behavior of K (e) differs slightly from the one- 
dimensional e~ 3 / 2 power law, but this deviation 
goes in the wrong direction to explain the discrep- 
ancy between theory and experiment. 

4.3. Diffusive approximation in narrow wires 

The energy scales probed by WL and Relax ex- 
periments are rather different. In wires, the value 
of t v is essentially determined by the low energy 
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Fig. 5. Energy dependence of the kernel 
K(e) of Coulomb interaction in a wire with 
L = 10 fim,w = 130 nm,t = 45 nm and D = 200 cm 2 /s. 
The asymptotic zero-, one- and three-dimensional regimes 
(0D, ID, 3D) arc characterized by K(e) = K(0), 
K(s) oc e~ 3 / 2 and K(e) oc e -1 / 2 , respectively (straight 
lines). The two-dimensional regime is not clearly vis- 
ible because w t. The range of relevant e's for 
the Relax experiments is determined by fcgT m ; n and 
et/max- The normalization factor on the y-axis is 
K(0) = {AbTvihD/L^rwpwtLy 1 



cut-off of the interaction, at h/r v . In the samples 
presented here, t v ranges (in the relevant temper- 
ature range: 1 K down to 40 mK) from 1 to 20 ns, 
corresponding to energies K/t v between 0.03 and 
0.6 jieV. In the Relax experiments, the shape of 
f(E) is entirely determined by energy exchanges 
of an amount between fc/jT and ell, in practice be- 
tween 4 and 500 fieV. According to Eq. (6), the 
characteristic lengthscale 1/q = \J hD/e for the in- 
teraction is therefore a few micrometers for WL, 
several hundreds of nanometers for Relax. The dis- 
crepancy between the results of the two types of 
experiment could point out a failure of the diffusive 
model, in which the QP dynamics is described by 
a single diffusion constant D. This argument is re- 
inforced by the fact that the elastic mean free path 
deduced from D is of the order of the wire thick- 
ness t, indicating that surface and grain boundary 
scattering dominate the elastic processes. If surface 
scattering alone were dominant, the elastic mean 
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free path of QPs with an instantaneous wavevec- 
tor along the axis of the wire would be very differ- 
ent from that of QPs travelling in a perpendicular 
direction, and the diffusive approximation would 
break down. To our knowledge, Coulomb interac- 
tion has never been investigated in this regime. 
However it is not clear why this situation could be 
described by the same energy dependence and why 
the intensity could be larger. 

4.4. Departure from equilibrium 

WL experiments are performed very close to 
equilibrium. In Relax experiments, a voltage U ^> 
ksT/e is applied to the wires in order to establish 
an out-of-equilibrium situation. Near the Fermi 
level, the distribution function is very different 
from a Fermi function, and it could be argued that 
the derivation leading to the expression (4) of the 
prefactor n ee is no longer valid. In order to test this 
hypothesis, we have performed a complementary 
experiment, described below, in which the effect 
of the distance to equilibrium is investigated. 



5. A new Relax experiment close to 
equilibrium 




Fig. 6. Schematic diagram of an experiment to measure 
fx{E) in a wire close to equilibrium. Quasiparticles are 
injected into the wire from a superconducting wire (labelled 
injector) through a tunnel junction biased at potential U 
(switch position 2). The distribution function f x (E) at 
position x = 0.25 is then determined from the dl / dV 
characteristic of the probe junction. Alternatively, the wire 
can be driven far from equilibrium by applying the voltage 
bias U across the wire (switch position 1). The resistance 
Rg is chosen so that the potential of the right reservoir 
remains close to zero when the switch is in position 2. 




E (meV) 

Fig. 7. Measured (o) distribution function fi (E) in the 

4 

"conventional" Relax experiment using sample AgXII40 
with the switch of Fig. 6 in position 1, and for U = 0.2 mV. 
The solid line is a numerical solution to the Boltzmann 
equation using the prefactor K ee = 0.18 ns — ^eV - 1 / 2 for 
the Coulomb interaction between electrons. As shown by 
the three dot-dashed lines, other values of ft ee produce 
markedly worse fits to the data. In particular, the theoret- 
ical value K ee = 0.025 ns~ 1 meV~ 1 / 2 does not come close 
to reproducing the experimental results. 

Fig. 6 shows a schematic of sample AgXII40, 
which was designed to investigate the effect of the 
deviation of f(E) from an equilibrium Fermi dis- 
tribution. As in other Relax experiments, a wire 
(38 /im long, 180 nm wide, 45 nm thick) is placed 
between large contact pads. A superconducting 
probe electrode is placed at x = 1/4, with a tun- 
nel resistance to the wire of 15 fcfi. The size of the 
tunnel junction was 0.18x0.23 /jm 2 . When the 
switch on Fig. 6 is placed in position 1, the "con- 
ventional" Relax experiment can be performed. A 
measured distribution function is shown in Fig. 7. 
The intensity of the Coulomb interaction deduced 
from the fits of f(E) is n ee = 0.18 ns^meV- 1 / 2 , 
as indicated in Table 4. Eq. (4) has been used 
[26] to calculate the theoretical value K*e y — 
0.025 ns -1 meV -1 / 2 . This discrepancy is of the 
same type as the one observed in the other sam- 
ples of Table 4. A second superconducting elec- 
trode, denoted injector in Fig. 6, forms a tunnel 
junction with the wire around its center, but with 
a much smaller resistance i?i n j = 1.1 kf2 than 
the probe junction, resulting from a larger area: 
0.57x0.8 £im 2 . This junction was obtained at the 
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overlap between the u>; n j = 0.8 /im-wide super- 
conducting electrode and the wire, which presents 
an intentional broadening at this position. When 
the switch of Fig. 6 is placed in position 2, quasi- 
particles are injected through the tunnel junction 
into the wire when \U\ > A/e, with A the gap in 
the QP density of states of the injector. On the 
normal side of the tunnel junction, the QP distri- 
bution function is therefore expected to display 
a step, the shape of which reflects the BCS den- 
sity of states n s {E) = Re (|£;| /VE 2 - A 2 ). The 
height of the step away from the BCS peak is given 
by the ratio of the injection rate of QPs to the 
diffusion rate towards the two normal reservoirs: 
fi(E) ~ (f ) /i? inj = r (the factor 1/4 results 
from the parallel combination of the two halves of 
the normal wire as will be shown below) . A quan- 
titative description follows from the introduction 
of new boundary conditions in the Boltzmann 
Eq. (3): f x (E) is a Fermi function with a zero 
electrochemical potential at x — and —eU r at 
x = 1, whereas at x — \ current conservation at 
each energy implies 



vpwteD 



with 



df x (E) 
Ldx 



df x {E) 
Ldx 1 



eRh 



-n s (E + eU)(f s (E + eU) - ft(E)) 



where fs(E) is the distribution function in the su- 
perconducting injector. We neglect here the slight 
modification of the DOS in the wire due to prox- 
imity effect, because of the small transparency of 
the tunnel barrier. Finally, 



df x (E) 

dx 1 
R 



df x {E) 
dx 



Ra 



n s (E + eU)(f s (E + eU) - h(E)). (7) 



The electrical potential of the right reservoir, which 
is connected to ground by a bias resistance Rb = 



12 SI, is given by U r 



1 rr b 



U. Since 



mation U r - 
U(E) = h- 



2Ri„ 



2 R+R B J " m J\ 

0.005, we make the approxi- 



0, so that the situation is symmetric: 
X (E) and Eq. (7) becomes 



df x (E) l 



df x {E). 



dx x ~i dx x ~? 

= 2rn s (E + eU)(f s (E + eU) - fi(E)). (8) 

In the absence of interactions, at T = 0, one ob- 
tains directly for x < \ (assuming U < —A): 

1 for E < 

f x (E) = <( 2xfi {E) for E € [0, -eU - A] 
for E > -eU - A 



and 



r n s {E + eU) 



(9) 



1 + rns (E + eU)' 
The spatial dependence of f x (E) is plotted in 
Fig. 8 for x < i, assuming r — 0.1 for visibility (in 
the experiment, r ~ 0.025). It is seen that f x (E) 
is much closer to a Fermi function than when the 
voltage is applied across the wire. 




Fig. 8. Schematic diagram showing the spatial and en- 
ergy dependence of the distribution function f x (E) of QPs 
driven out-of-equilibrium by the voltage U using the ge- 
ometry with the switch of Fig. 6 in position 2 (we have as- 
sumed U < —A). The surrounding box shows the density 
of states along the circuit and the gray volume shows the 
occupied states whose normalized density is f x (E)n x (E). 
The inelastic processes involving QPs are assumed to be 
very weak for clarity. The thick line shows the distribution 
function fi (E) at x = 1/4. 

4 

An experimental curve, obtained for U — 
—0.27 mV, is shown in Fig. 9. As predicted, it 
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f(E) 




-0.05 0.00 0.05 0.10 0.15 



E (meV) 

Fig. 9. Measured (o) distribution function fi (E) in the new 

4 

Relax experiment using sample AgXII40 depicted in Fig. 6, 
with the switch in position 2, and for U = —0.27 mV. 
The data are also shown magnified by a factor 10 (right 
scale). The solid line is a numerical solution to the 
Boltzmann equation with boundary condition given by 
Eq. (9), using as prefactor for the Coulomb interaction 
K ee = 0.18 ns~ 1 meV~ 1 / 2 . The two dot-dashed lines show 
that other values of re ee produce markedly worse fits to the 
data. 



presents a very small step (fi(E) w 0.025) ex- 
tending from E = to E = —ell — A, with 
A = 0.18 mV the gap for the injector deduced 
from its I — V characteristic, measured separately. 
The blow-up (x 10, right scale) shows the expected 
small peak near E = —ell— A. We also show f(E) 
calculated using the same parameters as those 
deduced from the "conventional" measurement, 
using Eq. (3) and (8). Except for a slight rounding 
of the small peak, the agreement is within experi- 
mental accuracy for all the values of U for which 
data were taken (—0.22 to —0.31 mV). We show in 
particular that other values of n ee would produce 
curves which significantly differ from the measured 
one. Hence the value of n ee deduced from energy 
exchange experiments does not seem to depend on 
whether the distribution is far from equilibrium, 
as in the original experiment (Fig. 7), or close to 
equilibrium, as in the newer experiment described 
here. Our conclusion is that Coulomb interaction 
is not modified by the fact that f(E) is not exactly 
a Fermi function. 



6. Conclusions 

In Section 4.1, we discussed the possibility that 
the anomalously high rates of energy exchange ob- 
served in many Relax experiments could be caused 
by residual magnetic impurities. Two arguments 
against this hypothesis were: 1) it seems implau- 
sible that all samples used in Relax experiments 
contain impurities that are not present in any sam- 
ple used for localization experiments, since both 
kinds of samples were fabricated in the same ap- 
paratus; and 2) we checked whether adding a term 
of the form K(e) oc e~ 2 to the interaction ker- 
nel could significantly decrease the value of n ee 
obtained from fitting the data to the solution of 
Eq. (3). But those two arguments do not rule out 
another possibility, namely that both kinds of sam- 
ples contain magnetic impurities with integer spin 
and with a magnetic anisotropy of the form JCS^ 
in the impurity Hamiltonian [29]. Such a term is 
predicted in the presence of spin-orbit scattering, 
for magnetic impurities located close to the sample 
surface [30] . If the characteristic energy JC satisfies 
ksT <C K < eU, then such impurities would con- 
tribute to energy exchange but not to dephasing. 
The contribution to K(e) from such impurities de- 
pends on both K, and B, but is not expected to 
be of the form K(e) oc e~ 2 . In principle, the pres- 
ence of such impurities should be detectable in ex- 
periments in the presence of a magnetic field. In- 
deed once g[iB 3> eU, their contribution vanishes. 
The absence of visible magnetic field dependence 
in sample AgX20 seems to rule out this possibility. 

In conclusion, the energy dependence of 
Coulomb interaction in disordered wires is well ex- 
plained by theory. The intensity of the interaction, 
as deduced from phase coherence time measure- 
ments, is quantitatively in agreement with theory, 
whereas for energy relaxation, an unexplained dis- 
crepancy remains. A new version of the Relax ex- 
periment has demonstrated that this discrepancy 
is not due to the out-of-equilibrium situation. 
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